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ABSTRACT 

We present proper motions for 21 bright main shell and 17 faint, higher-velocity, outer ejecta 
knots in the Cas A supernova remnant and use them to derive new estimates for the remnant's 
expansion center and age. Our study included 1951 — 1976 Palomar 5 m prime focus plates, 
1988 - 1999 CCD images from the KPNO 4 m and MDM 2.4 m telescopes, and 1999 HST 
WFPC2 images. Measurable positions covered a 23 to 41 yr time span for most knots, with a 
few outer knots followed for almost 48 yr. We derive an expansion center of q;(J2000) — 23*^ 
23™ 2T.n ±0!05, (5(J2000) = 58° 48' 49'.'4 ±0'.'4 (ICRS), with little difference between centers 
derived using outer or main shell knots. This position is 3'.'0 due north of that estimated by 
van den Bergh & Kamper (1983). It also hes 6'.'6 ±1'.'5 almost due north (PA = 354°) of the 
remnant's recently-detected central X-ray point source, implying a transverse velocity for the 
X-ray point source ~ 330 km s^^ at a distance of 3.4 kpc. Using the knots which lie out ahead 
of the remnant's forward blast wave, we estimate a knot convergent date of A.D. 1671.3 ±0.9 
assuming no deceleration. However, a deceleration of just ~ 1.6 km s~^ yr""'^ over a 300 yr time 
span would produce an explosion date ~ A.D. 1680, consistent with the suspected sighting of the 
Cas A supernova by J. Flamsteed. 

Subject headings: ISM: individual (Cassiopeia A) - supernova remnants - ISM: kinematics and dynamics 



1. Introduction 

Cassiopeia A (Cas A) is the youngest Galactic 
supernova remnant (SNR) known and, with the 
exception of the Sun, ranks as the strongest dis- 
crete radio source in the sky at 100 — 1000 MHz. 



^ Based in part on observations with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc. under NASA 
contract No. NAS5-26555. 



At radio, optical, and X-ray wavelengths, Cas A 
consists of a c± 2' radius broken shell of SN debris 
rich in O, S, Si, Ar, Ca expanding at 4000 - 6000 
km s~^. Within this shell lie about two dozen 
knots of much slower moving, N-rich clumps of 
pre-SN, circumstellar mass loss material. Outside 
of the shell, faint radio and X-ray emission ex- 
tends to a radial distance of ~ 160" where a faint, 
filamentary edge of X-ray emission marks the cur- 
rent location of the remnant's forward shock front. 
At an estimated distance of 3.4j^Q j kpc (Reed et 
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al. 1995), these angular dimensions correspond to 
main shell and outer shock front radii of 2 pc and 
2.7 pc respectively. Several dozen faint optical 
knots with velocities of 8000 to 15,000 km have 
been detected outside some sections of the main 
shell, mainly in a northeastern "jet" of high speed 
ejecta (Fesen & Gunderson 1996; Fesen 2001). 

The remnant's precise age is uncertain. From 
proper motion studies for ~ 100 of Gas A's op- 
tical knots during 1951 1980, Kamper & van 
den Bcrgh (1976, hereafter KvdB76) and van den 
Bcrgh & Kamper (1983, hereafter vdBK83) deter- 
mined an explosion date of 1658 ±3 for the rem- 
nant as a whole (assuming no deceleration) and 
a somewhat later date of 1671 for a few higher- 
velocity northeastern "jet" knots. The difference 
between these derived dates probably reflects a 
greater deceleration of bright main shell knots 
caused by their interaction with the remnant's re- 
verse shock. 

There are no unambiguous historical observa- 
tions of a bright nova or variable star in Gas- 
siopeia that might be associated with a late 17-th 
century supernova. However, on 1680 August 16 
John Flamsteed, the first Astronomer Royal, re- 
ported seeing a 5th 6th magnitude star he desig- 
nated "supra r" and later renamed 3 Cassiopeiae 
in his 1725 Historia coelestia star catalog (Ash- 
worth 1980). Its proximity to Gas A, together with 
the fact that he never observed this star again, 
raises the possibility that he sighted the Gas A 
supernova in the summer of 1680. 

The positional differences between Gas A and 
the 3 Gas position are, however, troublingly large. 
Flamstced's location for 3 Gas is offset from Gas 
A by 12.'1 in right ascension and 8'.6 in declina- 
tion. Although refraction and sextant corrections 
might decrease these residuals to ~ 6' in both coor- 
dinates (errors not unprecedented for Flamsteed), 
the case for Flamsteed's sighting of Gas A is con- 
troversial (Broughton 1979; Kamper 1980; Hughes 
1980). 

Without additional evidence, the significance of 
Flamsteed's observation might well remain incon- 
clusive. However, the large proper motions of the 
remnant's ejecta knots (/x = 0'.'4 — 0'.'6 yr~^) can 
be used to set limits on Gas A's age and thereby 
test the possible 1680 explosion date. An accu- 
rate measurement of Gas A's age would in turn 
provide key information about deceleration of its 



high-speed knots and thus the phase of its evolu- 
tionary development. 

Accurate proper motion measurements can also 
be used to improve determinations of the the rem- 
nant's center of expansion. KvdB76 determined 
the remnant's expansion center to within an error 
radius of about one arcsecond. No optical point 
source or extended emission is present at this lo- 
cation down to I ~ 24 mag (van den Bergh & 
Pritchet 1986). 

Knowledge of Gas A's precise expansion cen- 
ter has recently gained greater interest with the 
Chandra X-Ray Observatory discovery of an X-ray 
point source near the remnant's center (Tanan- 
baum 1999). This object, which could be either a 
neutron star with magnetized polar caps or an ac- 
creting compact object (Pavlov et al. 2000; Umeda 
et al. 2000; Ghakrabarty et al. 2001), lies sig- 
nificantly offset from estimates for the remnant's 
center of expansion (GOE). Point source coor- 
dinates derived from ROSAT and Chandra data 
show a separation of 1" — 5" from the KvdB76 and 
vdBK83 GOE, and some 16" - 20" from the GOE 
inferred by Reed et al. (1995) using knot radial ve- 
locities. These offsets imply transverse velocities 
of 50 - 250 km s'^ and 800 - 1000 km s"! respec- 
tively, assuming d = 3.4 kpc and an age of 320 yr 
(Pavlov et al. 2000). 

In this paper, we present proper motions of 17 
outlying high- velocity ejecta knots discovered over 
the last decade [see Fesen (2001) and references 
therein] along with 21 selected main shell knots. 
Many of these 40 knots can be seen on the earli- 
est archival Palomar 5 m PF plates, giving proper 
motion baselines of nearly five decades, or about 
l/7th of the remnant's age. We use these proper 
motions to determine a more accurate position and 
date for the supernova. 

2. Ejecta Knot Observations 
2.1. Image Data 

Our observational material includes Palomar 
5 m prime focus (PF) plates dating back to 1951, 
GGD images taken 1988 - 1999 with the KPNO 
4 m and MDM 2.4 m telescopes, and a few 1999 
EST WFPG2 images. Table 1 lists information 
on the images that were used. Although archival 
plates and modern GGD image data have, in many 
cases, substantially different spectral sensitivities, 



2 



the angular scale of knot emission stratification 
or ionization structures lies well below all but the 
highest resolution HST data (Fesen et al. 2002) 
and thus docs not pose a significant problem for 
the inter-comparison of knot positions from these 
different data sets. 

To maximize the time base, wo examined sev- 
eral dozen archival Palomar 5 m PF plates begin- 
ning with R. Minkowski in 1951 and ending with S. 
van den Bcrgli in 1989. Most of these plates were 
unsuitable for this project because of poor image 
quality or weak knot detection, but four were se- 
lected for use. These included two plates taken on 
back-to-back nights in 1951, a better and deeper 
1958 image (van den Bergh & Dodd 1970, their 
Fig. 1), and a superb 1976 image (van den Bergh 
& Kamper 1983). 

Modern CCD interference and broadband filter 
images of all or portions of the Cas A remnant 
obtained from November 1988 through October 
1999 were also measured. Several of these have 
been used in prior studies, and a few were taken 
expressly for this project. Some high resolution 
1999 epoch WFPC2 HST images were also used 
for several outlying northern and jet knots. 

We selected for measurement a total of 38 ejecta 
knots that had measurable positions on the Palo- 
mar 1976 and later images (At = 23 yr). Many 
were outer knots, together with some shell knots 
visible from 1958 through 1999 (At = 41 yr). 
Nearly a dozen knots, mostly among the outer 
ejecta, were detectable from 1951 through 1999, 
a span of almost 48 years, which is about 15% of 
Cas A's age. 

It is unlikely that suitable earlier images exist. 
The Palomar 5 m was completed in 1948, coinci- 
dent with the discovery of Cas A as a localized 
radio source (Ryle & Smith 1948). Furthermore, 
the remnant has brightened significantly over the 
last half-century (van den Bergh & Kamper 1985), 
and the individual knots often have finite visibility 
lifetimes (KvdB76). These factors set a practical 
limit of ~ 50 yr for the time span over which the 
proper motion can be studied at this time. 

2.2. Knot Selection 

Of the 38 knots selected, 21 were main shell 

features and 17 were outer, higher-velocity knots 
in the NE jet or elsewhere. Figure 1 shows the 



Fig. 1.— [S II] AA6716,6731 image of Cas A from 
1992 July (Fesen & Gunderson 1996). The lower 
panel shows all the knots used in this study. The 
upper panel is a magnified view of the main shell 
region. The positions marked are derived from 
the knot trajectories adjusted to the epoch of the 
image. Some of the knots are not visible. 

location of the selected knots. Tabic 2 cross-lists 
our designations with those of earlier studies where 
possible. 

Our two main knot selection criteria were: (1) 

Distinct appearance with an absence of signifi- 
cant morphological changes. We looked for knots 
which were compact with steady morphologies 
that allowed secure identification and centroiding. 
None of the knots was perfect but some, like the 
bright outer Knot 15, provided excellent positional 
measurements over the entire 48 yr time span 
surveyed. Like most other outlying fast-moving 
knots, it has a relatively stable morphology (see 
Fesen 2001). On the other hand, main shell knots 
can show substantial changes in appearance on im- 
ages separated by just 5 to 10 yr. Consequently, 
we included only the most distinct and persistent 
main-shell features. The number of main-shell fea- 
tures we used was therefore relatively small com- 
pared to prior studies (KvdB76 and vdBK83). (2) 
Time span of visibility. The longer a knot is mea- 
surable, the greater weight it has for determining 
the remnant's expansion center. We therefore bi- 
ased our knot selection toward knots with long 
visibility time spans (> 20 yr). This resulted in a 
much smaller knot sample than the 102 measured 
by KvdB76. Their knots covered time spans rang- 
ing from 3 — 24 yr, with 46% of their knots visible 
for less than 15 yr. 

3. Astrometric Procedures and Measure- 
ments 

3.1. Reference Star Grid 

We began by constructing a grid of reference 
stars. First, centroids for several hundred unsat- 
urated stars were measured on the Ha and [S II] 
AA6716,6731 images taken with the MDM 1.3 m 
telescope in 1992 July, using the IRAF incarna- 
tion of DAOPHOT (Stetson 1987). We then cross- 
identified these stars with the USNO A2.0 cata- 
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loguc (Monet ct al. 1996) and derived a 6-constant 
linear plate model, rejecting stars with large resid- 
uals. Because the centroids in the CCD data have 
much better internal precision than the USNO 
A2.0 coordinates, we transformed the CCD cen- 
troids to right ascension and declination using the 
plate model, and averaged the results from the 
Ha and [S II] exposures. This procedure yielded 
right ascensions and declinations approximately 
on the International Coordinate Reference System 
(ICRS) of the USNO A2.0, but with much higher 
internal precision. The USNO A2.0 is based on the 
original Palomar Observatory Sky Survey plates, 
which for this field are epoch 1954.6, so the elim- 
ination of high-residual stars removes stars with 
appreciable proper motions, as well as blended im- 
ages and other difficult cases. For the final step, 
the refined celestial coordinates were converted to 
tangent-plane coordinates, using the KvdB76 cen- 
ter of expansion as the tangent point. The final 
reference grid consisted of 141 stars with red mag- 
nitudes (from the USNO A2.0) from 15.2 to 19.1 
covering an 8' x 8' field. 

Because this reference star grid is fundamental 
to all our results, we checked it for various sources 
of error: (1) errors in the grid's coordinate zero 
point which would affect comparisons with other 
results, (2) proper motion of the grid stars, and (3) 
radial distortions in the grid which could cause a 
systematic error in the age estimate. Below, we 
consider each of these sources of error in turn. 

3.1.1. Zero Point for Reference Grid 

We performed two checks on the positional zero 
point of the reference star grid. 

First, we examined the twelve reference stars 
(for epoch 1965) tabulated by KvdB76. We mea- 
sured these stars on all the images on which they 
appeared, used our reference star net to derive po- 
sitions and proper motion as described below for 
the knots, and compared the results to those tabu- 
lated in KvdB76. Ten of the 12 KvdB76 reference 
stars had sufficient observations in our data. For 
these, our right ascensions were on average 485 
mas larger than theirs, with an RMS scatter of 67 
mas, while our declinations were 448 mas smaller, 
with a scatter of 112 mas. 

Second, during the preparation of this work, 
the Tycho-2 astrometric catalogue became avail- 



able (H0g et al. 2000). Four Tycho-2 stars ap- 
pear on enough images to derive good positions 
for epoch J2000, and for three stars the time base 
gives adequate proper motions. Because the stars 
were highly saturated in nearly all our pictures, we 
estimated the star' centers and their uncertainties 
by eye. Our derived epoch 2000 positions for the 
four Tycho stars are, on average 121 ± 80 mas east 
and 82 ± 79 mas south of the catalogue positions, 
well within the accuracy to which the USNO A2.0 
is expected to align with the ICRS. 

The Tycho and USNO catalogues are based 
on Hipparcos observations, which should provide 
much more reliable all-sky positions than the cat- 
alogs available to KvdB76. We conclude that our 
grid is in registration with the ICRS to within 
±0'.'2 at worst, and that KvdB may have suffered 
a barely significant zero-point error. 

3.1.2. Proper Motions of Grid, Stars 

In deriving plate models, we did not adjust the 
positions of our grid stars for proper motions. We 
implicitly assumed that the proper motions of the 
grid stars were small. Our grid stars are on aver- 
age only a little brighter than the dozen y ~ 19 
mag reference stars used by KvdB 76. They re- 
mark that stars this faint should have intrinsic 
proper motions from the solar motion and Galac- 
tic rotation of ^ 1 mas yr"-'^. It is therefore likely 
that our grid stars' motions are similarly small. A 
systematic offset Hsys in proper motions can seri- 
ously affect the derived center for the SNR, since 
it displaces the center by /UgysX ~ 300 yr. 

For the three Tycho-2 stars for which proper 
motions could be derived on our reference grid, 
the weighted averages of /xt2 — Mgrid were —3.5 ± 
2.5 mas yr~^ in right ascension and -1-0.8 ± 2.5 
mas yr~^ in declination. Therefore, there was no 
evidence for a significant motion with respect to 
the ICRS. 

The comparison with the KvdB reference stars 

was a bit more complex. They explicitly derived 
proper motions for their reference stars, and found 
a 9 mas yr~^ mean proper motion in declination. 
They decided this was spurious and later adjusted 
their derived Cas A center to account for the drift 
in their reference grid. Reducing 10 of their ref- 
erence stars with respect to our grid gives mean 
differences /zkv — A*ref of —2.8 ± 1.4 mas yr~^ in 
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right ascension and +8.3 ±1.4 mas yr"-'^ in dec- 
lination. This is just as expected on the basis of 
their remarks. Indeed, our procedure is logically 
similar to theirs, but less complex: we simply as- 
sumed the faint grid to be motionless ab initio, 
while they derived proper motions and corrected 
them at the end to achieve the same result. 

Finally, we checked our grid stars individually 
against the reference grid, deriving proper motions 
for each star based on all the images on which it 
appeared. Formally, this was flawed by the inclu- 
sion of the star itself in the plate models, but with 
our sample of 140 stars, the effect should be neg- 
ligible. Grid stars observed over the whole range 
of dates had typical estimated proper motion un- 
certainties of ± ~ 1.7 mas yr""*-, and relatively 
few showed significant proper motions. In the 
fitting procedure used for the images (described 
below), grid stars with largo residuals were itera- 
tively clipped out, so the few stars with significant 
proper motions did not affect our results. 

Prom all these tests, we conservatively estimate 
the grid to be inertial to within 2 mas yr~^. 

3.1.3. Field Distortions 

The procedure used to set up the reference grid 
is valid provided that distortions in the field of 
the MDM 1.3 m telescope are insignificant. This 
is likely to be the case. Cudworth & Rees (1991) 
measured the field distortions of several southern 
telescopes, including the CTIO 1.5 m which, like 
the MDM 1.3 m, is an f/7.5 Ritchcy- Chretien re- 
fiector. If one equipped the CTIO 1.5 m with a 
CCD having the same size as that used in deriv- 
ing our reference grid, then their radial distortion 
term ag would contribute only 7 mas at the corners 
of the field of view. 

Nonetheless, we searched for field distortions in 

several different ways. (1) In fitting the USNO 
A2.0 stars, we did not see any trends in the resid- 
uals from the 6-constant plate model. The USNO 
A2.0 typically has centroiding errors in the 250 - 
500 mas range. In view of the number of stars 
used, this alone limits systematic trends to < 200 
mas. (2) We examined archival CCD images of 
the globular cluster M13 taken with the MDM 1.3 
m and the same camera as the Cas A images. Dr. 
Kyle Cudworth kindly provided us with a list of 
star positions in M13, which he estimated were 



accurate to ~ 20 mas for relative positions. A 
six-constant fit of the CCD centroids to those posi- 
tions gave an RMS residual of 70 mas, again show- 
ing no obvious systematic trends. When a more 
elaborate model was used, the residuals were not 
improved significantly. The scatter is somewhat 
larger than expected, but docs not seem to indi- 
cate any field distortions. (3) In 1999 October we 
obtained a set of short 7-band CCD exposures of 
Cas A with the MDM 2.4 m telescope, covering an 
8' X 8' field. A 6-constant fit to the 124 standard 
reference stars included in this image gave an RMS 
residual of 60 mas, without any iteration; iterative 
clipping of high residuals brought this down to 34 
mas, with 110 stars remaining. Again, a more 
elaborate fit did not result in significant improve- 
ment. 

Because the reference grid is based on expo- 
sures taken with the MDM 1.3 m, this last test 
simply compares the two telescopes. However, in 
January 2000 we also performed an astrometric 
calibration of the 2.4 m by obtaining two sets of 
short-exposure 2.4 m images of a portion of the 
Stone et al. (1999) astrometric standard region E. 
In one set of images, we fitted 172 stars to a 6- 
constant plate model. This gave a 153 mas RMS 
residual. Iterative elimination of the largest resid- 
uals brought this down to 53 mas with 114 stars. 
Similar results were found with the other set of 
images. The residual maps of the two sets of im- 
ages did not show systematic distortions, but were 
highly correlated with each other. This suggests 
that most of the error arises from the catalogue po- 
sitions, probably due to the (necessary) inclusion 
of many stars near the faint limit of the catalogue. 

In summary, the tests we made did not show 
any geometric distortion in our reference grid. 
The results suggest that systematic distortions are 
smaller than ~ 60 mas and that the centering pre- 
cision of the reference stars is conservatively ~ 50 
mas. 

3.2. Image Solutions 

We scanned the four Palomar PF plates (Ta- 
ble 1) on the Yale Astronomy Department's PDS 
microdensitometer. We used a 13.3 fj,m x 13.3 fim 
scanning aperture and sampled every 12.656 fxm in 
a 3300 X 4100 raster centered on the remnant, the 
long dimension being east- west. The plate scale 
was 11'.' 1 mm~^, yielding 0'.'141 pixel" ^. We used 
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SExtractor (Bertin & Arnouts 1996) to derive star 
centers on for the Palomar plates. 

The Palomar 5 m prime-focus camera used a 
corrector which produced substantial radial dis- 
tortions. Similar corrector distortions have been 
discussed by Murray (1971) and Cudworth fc Rees 
(1991). Plate coordinates of the distortion cen- 
ter, (xo,t/o), enter the least-squares models in a 
non-linear fashion, so following Murray (1971) and 
Cudworth & Rees (1991) wc estimated (xo.yo) in 
a separate step. We first fit the reference stars 
with a simple 6-constant plate model of the form 

Xo = ao + aix + a2y, 

where Xq is the standard reference star coordi- 
nate, and X and y are the coordinates on the Palo- 
mar plate. This fit gave root-mean-square (rms) 
residuals of ^ 240 mas, with maximum values 
~ 690 mas. The distortion center {xo,yo) could 
then easily be estimated from the residual maps. 
These were then used in a 16-constant model of 
the form 

^0 = ao+aix+a2y+a3X^+a4xy+a5y'^+aexr'^+aTxr* , 

where = ((.x — a;o)^ + (y — j/o)^)"'^^^. This model is 
similar to that used by Cudworth & Rees (1991), 
but without magnitude terms. As expected, these 
fits were much better; after a few of the highest- 
residual stars were rejected, rms residuals ranged 
from 71 mas (P7252) to 133 mas (P553B). 

Fitting the CCD images was more straight- 
forward. We centroided the reference stars with 
DAOFIND and matched them to their standard 
star XY coordinates. For the ground-based CCD 
images, there were always > 20 stars matched. 
Because distortions in the CCD images were ex- 
pected to be relatively small, we used the 6- 
constant model described above, with any unruly 
stars omitted by iterative clipping. 

Reduction of the few HST images was not as 
simple. First, the images from the four CCDs 
of the WFPC2 were interpolated onto a single 
grid using the imnosaic task in the IRAF STS- 
DAS package, which approximately corrects for 
field distortions. Unfortunately, the WFPC2 field 
of view is so small that one of our fields had only 
six reference stars. To increase the number of ref- 
erence stars we measured some fainter stars on 
the 'wide [S II]' 2.4 m image from 1999 October, 



transformed these over to the standard grid and 
used them in the fit to the HST data. With these 
added stars, each HST field had at least 14 ref- 
erence stars. The rms scatter for the 6-constant 
models was 50 to 55 mas. Holtzman et al. (1995) 
quote an rms scatter of 10 mas in their fit to the 
WFPC2 field distortions. The larger scatter found 
here probably arises from the errors in the refer- 
ence star positions. In particular, there was no 
pattern indicating that the zero-point offsets be- 
tween the CCDs were different from those assumed 
by wmosaic. 

3.3. Knot Measurements and Fits 

We tried several methods to measure positions 
of the selected 38 knots. This was a somewhat 
complicated problem, since many of the knots were 
resolved in our images. The dominant source of 
centering uncertainty resulted from structure in 
the knots, rather than photon or grain noise. In 
the end, we simply estimated the knot centers by 
eye using a cursor on an image display. We al- 
lowed our judgment to be informed by centroids 
from such tools as the imexamine task in IRAF 
and (for the photographic scans) the centers from 
sextractor. The selected knots often had 'head- 
tail' structures, in which case we estimated the 
center of the head. Some of the shell knots were 
embedded in nebulosity. In such cases, we tried to 
center on the brightest part of the knot. 

Our uncertainty estimates were also subjec- 
tive, but we attempted to err on the conservative 
side. When we later fit straight-line trajectories to 
the knots, we found that the residuals were often 
smaller than one would expect on the basis of our 
estimated uncertainties, demonstrating that our 
estimates were indeed conservative. For the photo- 
graphic images, our estimated errors were guided 
by lower bounds on the error based on the noise in 
the plate fog, the central brightness of the knot, 
and the knot's angular size. These bounds were 
calibrated with a Monte Carlo simulation. 

Once xy coordinates were measured on all the 
images, they were transformed to the tangent- 
plane coordinate system using the solutions de- 
scribed earlier. The rms scatter in the plate solu- 
tion was added in quadrature to each the knot's 
estimated uncertainty. If the rms scatter was less 
than 40 mas, it was set to 40 mas to account for 
systematics. Transformed knot positions were au- 
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tomatically collated with the exposure epochs for 
the images, yielding a time series of positions for 
each knot. 

Finally, each knot's trajectory was fitted with 
a straight line, 

X{t)=Xo{to)+Hx{t-to), 

where to is the weighted mean epoch of observa- 
tion, and similarly for Y. Uncertainties were prop- 
agated into the coefficients in a standard manner. 
Table 3 gives the results of this procedure. 

4. Results and Discussion 

4.1. Cas A's Center of Expansion 

Figure 2a shows the trajectories of our selected 
knots extrapolated back to A.D. 1600, somewhat 
before the estimated explosion date. The dots arc 
the individual knot measurements with the width 
of each line indicating its statistical weight. 

Because the knots may have decelerated by 
varying amounts, we estimated the center using 
only the knots' lines of position, and did not ap- 
ply any constraints arising from the time; depen- 
dence (i.e, forcing the knots to start at the same 
epoch). We constructed a trajectory for each knot 
and computed its positional uncertainty near 
the time of the explosion by propagating the esti- 
mated position and proper motion errors. Because 
of the long time lever, the proper motion uncer- 
tainty dominated the errors in all cases. With this 
information we could compute, for any X and Y, 
a likelihood function of the form 

\{X,Y) = n ^ exp(-d2^/2<7?o), 

where dij_ is the perpendicular distance between 
{X, Y) and the knot's line of position. The {X, Y) 
which maximizes this is our estimate of the expan- 
sion center. 

This procedure gave X = +91 mas and Y = 
-1-2812 mas, referred to the KvdB76 center. This 
translates into a(J2GGG) 23'^ 23"" 27?77 ±0!05, 
(5(J2000) = 58° 48' 49'.'4 ±0'.'4. The center derived 
using only the outer knots nearly coincided with 
that derived from the selected shell knots. Table 4 
lists our main shell and outer knot centroids sep- 
arately along with our final values for the whole 
sample. The errors given are purely statistical. 



(a): Cas A — Knot Trajectories 
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(b): Cas A Knot Trajectories 
(Van den Bergli & Kamper 1976) 
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Fig. 2. — a) Trajectories of the knots used in this 
study. The dots are individual measurements, and 
the lines are fits to the trajectories. The widths 
of the lines increase according to the weight in the 
solution, b) Trajectories of the 102 knots used by 
Kamper & van den Bergh (1976). 



based on the Monte Carlo calibration of the cen- 
tering errors (see below). 

4.1.1. Error Estimates 

Because our expansion center differs from those 
of previous estimates (see below), we estimated 
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Fig. 3. Magnified view of the region in which 
the knot trajectories intersect. The straight hnes 
are the fitted trajectories, again with the line 
weight indicating tlic statistical weight. The con- 
tours represent confidence intervals for the cen- 
troid based on the likelihood ratio test, and repre- 
sent 95% to 99.5% confidence. The dots are cen- 
troids determined from Monte Carlo simulations of 
the measurement. A 'fudge factor' of 0.7 has been 
applied to the estimated errors for this calculation 
(see text). The origin of this plot is the KvdB76 
estimate of the center of expansion; note the sig- 
nificant offset of the new determination from this 
position. 

our measurement uncertainties in several ways. 

Since we have computed a likelihood function, 
the likelihood-ratio test described by Cash (1979) 
can be used to form confidence contours. Such 
contours, however, assume that the positional er- 
rors we estimated by eye arc truly one standard 
deviation. The resulting 95% confidence contour 
is an oval slightly elongated northwest-southeast, 
with a radius ^ 1'.'3. Figure 3 is a magnified view 
of the center region, showing the lines of position 
together with the 95% to 99.5% confidence con- 
tours from this procedure. 

We also used a Monte Carlo simulation of the 
proper motion measurements to check our uncer- 
tainties and to normalize our error estimates. We 
began by assuming that the explosion occurred at 



the observed maximum-likelihood position. Then, 
taking each knot's present position as known and 
fixed (because the proper motions dominate the 
errors), we computed an idealized proper motion 
for each knot, which extrapolated back exactly to 
the observed maximum-likelihood position. We 
next created 1000 artificial data sets by adding 
Gaussian random noise to the idealized proper mo- 
tions. For each artificial data set, we computed 
the maximum-likelihood center and its associated 
maximum A. 

In the first trials, the standard deviation used 
for the Gaussian noise was simply the estimated 
proper motion uncertainty of each knot - in effect, 
we took our estimated position errors to be realis- 
tic. For nearly all the trials, this choice led to max- 
imum values of A larger than observed, indicating 
that the errors in the proper motions were over- 
estimated (as expected since our error estimates 
were believed conservative). In the final Monte 
Carlo calculation, we multiplied the proper mo- 
tion standard deviations by a global 'fudge factor' 
/-0.7. 

These simulated data sets yielded maximum 
values of X{X,Y) very similar to that of the real 
data. The spread in the positions generated by 
this procedure should be a realistic indicator of the 
uncertainty in the centroid. The cloud of Monte 
Carlo positions is also shown in Figure 3 and has 
standard deviations in X and Y of 398 and 366 
mas, respectively. Half of the points lie within 
448 mas of the mean position and 95% within 946 
mas. This agrees well with the likelihood-ratio 
test described above once the scaling factor of 0.7 
is taken into accoTint and provides a cross-check 
on the likelihood-ratio procedure. 

As an additional check, we divided our data 
into two samples shell knots and outer knots. 
The outer knots gave X = 599 mas, Y = 2809 mas 
with half the Monte Carlo points within 508 mas of 
the mean; the shell knots gave X = —388 mas and 
Y = 2840 mas, with half within 783 mas. The dis- 
agreements between these and the combined data 
are about as expected given the estimated statis- 
tical errors. 

In a final statistical experiment we repeatedly 
selected half the knots at random found the max- 
imum likelihood center of the subsample. For this 
we used the observed trajectories without adding 
artificial noise. The distribution of the centers pro- 
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duced by 1000 such trials was rather elongated 
east- west {ax = 666 mas, ay = 410 mas), but did 
not extend toward a particular direction (e.g., the 
KvdB76 and vdBK83 center). This shows that 
our center determination is not thrown off by a 
few stray knots. 

4-1-2. Comparison to Previous Results 

Previous center of expansion values from op- 
tical and radio data are listed in Table 4. The 
first accurate estimate of Cas A's expansion point 
was made by van den Bergh & Dodd (1970) using 
Palomar 5 m PF plates of 27 "fast-moving knots" 
covering the time period 1951 through 1969. Kam- 
per & van den Bergh (1976) later updated this 
to include a total of 102 knots covering the addi- 
tional period 1970 — 1975. This study itself was 
supplemented by 1976 — 1980 measurements of 46 
especially long-lived knots by van den Bergh & 
Kamper (1983). The last two studies reached es- 
sentially the same central position estimate within 
measurement errors. More recently, Reed et al. 
(1995) found a significantly different expansion 
point based on a least-squares spherical fit to a 
plot of main shell knot radial velocities. However, 
this displaced center reflects radial velocity differ- 
ences between back and front hemispheres and is 
thus unlikely to be an accurate measure of the 
remnant's expansion center. 

Our estimated center lies only 3 arc seconds 
due north of the center derived by KvdB76 and 
vdBK83, but this is well outside their ±0'.'8 - I'.'O 
estimated uncertainty. Although our study in- 
cludes far fewer knots (38) than they used, our 
results have smaller formal errors and the trajec- 
tories have a smaller dispersion. For comparison 
we show in Figure 2b trajectories of the 102 knots 
measured in the KvdB76 study at the same scale 
as the trajectories in Figure 2a. The larger spread 
in knot trajectories in the KvdB76 data largely re- 
flects inclusion of fewer outlying knots and more 
main shell knots in their study. 

Finally, we note that the two reported centers 
derived from radio measurements lie significantly 
(~15") east of the centers estimated from the op- 
tical knots. Because the remnant's radio emission 
does not exhibit a smooth, globally-coherent ra- 
dial expansion (Bell 1977; Tuffs 1986; Anderson & 
Rudnick 1995), we believe that radio-derived cen- 
ters are less meaningful. 



Fig. 4. — Image of the center of the Cas A super- 
nova remnant taken with the MDM 2.4 m Hiltner 
telescope and a [S II] AA6716,6731 filter in 1" see- 
ing. The positions of the Chandra X-ray point 
source, the KvdB76 expansion center, and our re- 
vised expansion center are indicated. 

4-1.3. X-ray Point Source 

First-light Chandra observations of Cas A re- 
vealed the presence of a point-like X-ray source 
near the center (Tananbaum 1999). The source 
was subsequently confirmed through inspection of 
archival ROSAT (Aschenbach 1999) and Einstein 
data (Pavlov & Zavlin 1999). The X-ray source's 
position has been refined slightly using further 
C/iandra observations (Murray et al. 2001; Kaplan, 
Kulkarni, & Murray 2001). Table 5 fists this re- 
fined Chandra position, together with those esti- 
mated from investigations of ROSAT and Einstein 
HRI data. Figure 4 shows the Chandra position, 
and both KvdB76's and our expansion centers su- 
perposed on an optical image of the central region. 

The Chandra X-ray point source lies some 6'.'6 
south (position angle ~ 354°) of our derived center 
of expansion for Cas A. Our new expansion cen- 
ter actually lies farther away from the Chandra 
position than does the KvdB76 center. Both the 
ROSAT and iJmsfem observations, however, place 
the X-ray point source a few arc seconds farther 
to the west and north, thus closer to our expan- 
sion center. Assimiing a common origin for the 
ejecta knots and the point source, our center and 
the Chandra position together imply a transverse 
velocity of ~ 330 km at a distance of 3.4 kpc. 
This velocity is not unusual for a young pulsar, 
if this is indeed what it is (Umeda et al. 2000; 
Chakrabarty et al. 2001; McLaughlin et al. 2001). 

4.2. Explosion Date and Knot Decelera- 
tion 

Limits on the date of the Cas A SN explosion 
bear on two questions. First, could Flamsteed 
have seen it in A.D. 1680? Second, how much 
have the knots decelerated over the last 300 yr? 

The outer and shell knot trajectories give sig- 
nificantly different ages, in the sense expected if 
the shell knots have decelerated more than the 
outer knots due to their interaction with a re- 
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Cas A — Date Estimates 
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Fig. 5. — Dates at which the extrapolated 

knot trajectories pass nearest to the maximum- 
hkelihood center. The size of the symbol is in- 
versely proportional to the estimated uncertainty. 
Outer knots (squares) arc toward the top. and 
shell knots (circles) toward the bottom. Note the 
retardation and increased scatter of the shell. 

verse shock. This trend is clearly visible in Figure 
5, which shows the times at which the individual 
knots are computed to have passed closest to the 
center (assuming no deceleration); we will refer to 
this quantity as the crossing time. (Note: The cen- 
ter position used to compute the crossing times is 
derived from the whole data set. The knots move 
quickly enough that small variations in the place- 
ment of the center do not affect the trend.) 

Because of their smaller dciccilciration, the outer 
knots offer a better estimate of the date of the 
Cas A supernova than ejecta in the bright shell. 
A straight average of the 17 outer knots' cross- 
ing times yields an explosion date of 1671.3 ± 0.9, 
while the 21 main shell knots yield 1662 ± 1.7. On 
the face of it, the outer knot data indicate a date 
nine years earlier than Flamsteed's 1680 sighting 
of 3 Cas. However, a deceleration of only ~ 0.1 
mas yr~^ for the outer knots would change the 
date by ~ 10 years. In our best-observed knots 
the la uncertainty in the deceleration approaches 
this value, but none of the knots show significant 



deceleration. There is also no trend for the fitted 
acceleration vectors to be pointed inward toward 
the center. Thus we can neither directly detect 
nor disprove decelerations large enough to make 
Flamsteed's A.D. 1680 sighting coincide with the 
explosion. 

The dispersion in the knots' crossing times also 
affects the case for Flamsteed's 1680 sighting. As 
noted above, the mean crossing time is displaced 
from 1680. If in addition the dispersion in crossing 
times were small, then in order for Flamsteed to 
have seen the explosion the different knots would 
need to have suffered nearly identical decelera- 
tions. This would be unlikely given the typical in- 
homogeneity of the ISM. However, Figure 5 shows 
enough scatter in the crossing times that a Flam- 
steed sighting remains plausible. 

If we assume for a moment that the explosion 
date really was 1680 and that the knots have decel- 
erated uniformly ever since, we can then compute 
an implied deceleration for each knot. For the nine 
best shell knots (those with formal l-tr crossing 
time uncertainties less than 5 years) , these implied 
decelerations range from 0.04 mas yr~^ for Knot 
9 to 0.14 mas yr^^ for Knot 120, with a mean of 
0.10 mas yr"^. At 3.4 kpc, 0.10 mas yr"-'^ corre- 
sponds to a transverse acceleration of only 1.6 km 
s~^ yr~^ or a velocity change of some 2-5% over 
the age of the remnant. 

Detecting velocity changes ~ 1 — 2 km s~^ yr~^ 
in these faint, outer ejecta knots to test plausi- 
ble explosion dates would be difficult, but per- 
haps not impossible. The sudden brightening of 
Knot 19 along the remnant's western limb during 
the early 1970's (Fcscn 2001) suggests it may be 
decelerating significantly at the present epoch. A 
direct measurements of a knot's present-day decel- 
eration could, in principle, be used to explore the 
density of the interstellar or circumstellar medium 
around the remnant. But at present not enough is 
known about the knots' masses, dimensions, and 
structure to draw reliable conclusions about their 
environment from their decelerations. 

4.3. Explicitly Time-Dependent Estimates 

As a check on the center and explosion date es- 
timates, we also estimated these quantities jointly. 

Using the fitted knot trajectories, we stepped 
through a range of dates t around the explosion. 
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At each date, we computed the weighted mean po- 
sition of the knots, and then computed 



knots, implying a greater deceleration of the main 
shell. 



where di (t) is the angular distance of knot i from 
the mean center at date t. The estimated explo- 
sion date is then tmin, the date which minimizes 
S. The weighted mean position of the knots at 
tmin is an estimate of the explosion center. This 
procedure yielded t = 1671.5, X = +723 and 
Y = +2362 for the outer knots, t = 1659.8, 
X = +585 and Y = +1623 for the sheU knots, 
and t = 1669.8, X = +856, and Y = +2693 for 
the entire sample of knots. Because the knots in 
our sample are distributed non-uniformly around 
the remnant's periphery, differences in date trans- 
late into differences in position in a complicated 
way. The likely differential deceleration of the 
knots therefore makes this estimate less reliable 
than the previous estimate, which is based on lines 
of position. Nonetheless, the results are broadly 
similar - all the estimates put the best center sig- 
nificantly north of KvdB76's center, and the shell 
knots show evidence of deceleration. 

5. Conclusions 

Our new proper motions of 21 main-shell and 17 
higher-velocity, outer ejecta knots in Gas A leads 
to improved estimates of the center of expansion 
and the age. We find the expansion center to be 
a(J2000) = 23*^ 23" 27;77 ±0^05, (5(J2000) = 58° 
48' 49'.'4 ±0'.'4, with little difference between the 
centers derived using outer or main shell knots. 
This expansion point lies 6'.' 6 ±1'.'5 to the north 
of the recently recognized X-ray point source. If 
the point source originated in the explosion, the 
position offset implies a transverse velocity of ~ 
330 km s~^ at a distance of 3.4 kpc. 

Using the outer knots, most of which are in 
front of the main blast wave, we estimate a date of 
explosion of 1671.3 ±0.9 assuming no deceleration. 
However, interaction with local CSM/ISM should 
decelerate the knots. If the velocities have declined 
by only a few percent over the age of the remnant, 
the remnant age would be consistent with a sus- 
pected sighting of the supernova by J. Flamsteed 
in 1680. The age derived from the main shell knots 
is greater by 9 yr than that derived from the outer 
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Table 1 

Observational Material Used for Proper Motions 



Date 


Telescope 


Plate No./ 


Scale 


Emulsion/ 


Region 


Exposures 


(U.T.) 




Image ID 


( "/Pix) 


Bandpass 


Imaged 




(s) 


1951 Nov 01 


Palomar 5 m 


553B 


0.1406 


103aE + RG2 


Whole SNR 


1 X 


7200 


1951 Nov 02 


Palomar 5 m 


563B 


0.1406 


103aE + RG2 


Whole SNR 


1 X 


7200 


1958 Aug 11 


Palomar 5 m 


3033S 


0.1406 


103aF + RG2 


Whole SNR 


1 X 


5400 


1976 Jul 02 


Palomar 5 m 


7252vB 


0.1406 


098-04 + RG645 


Whole SNR 


1 X 


7200 


1988 Nov 10 


KPNO 4 m 




0.297 


[S II] A6725 


Jet, East Limb 


5 X 


128 


1988 Nov 10 


KPNO 4 m 




0.297 


Hq+[N II] 


Jet, East Limb 


5 X 


192 


1992 Jul 05 


MDM 1.3 m 




0.635 


broad [S II] A6725 


Whole SNR 


3 X 


1800 


1992 Jul 05 


MDM 1.3 m 




0.635 


Ha+[N II] 


Whole SNR 


3 X 


1800 


1996 Oct 06 


MDM 2.4 m 




0.275 


broad [S II] A6725 


NW.SW.NE.SE 


2 X 


1000 


1996 Oct 07 


MDM 2.4 m 




0.275 


Ha+[N II] 


NW,SW,NE,SE 


2 X 


600 


1996 Oct 07 


MDM 2.4 m 




0.275 


Cont. 6450 A 


NW,SW,NE 


2 X 


1000 


1999 Jun 12 


HST 2.4 m 


U52B0194R-0109R 


0.0996 


[S II] (F673N) 


Jet 


6 X 


1000 


1999 Jun 12 


HST 2.4 m 


U52B010CR-010FR 


0.0996 


R (F675W) 


Jet 


4 X 


600 


1999 Jun 13 


HST 2.4 m 


U52B0205R-0207R 


0.0996 


R (F675W) 


NW Limb 


3 X 


500 


1999 Jun 13 


HST 2.4 m 


U52B0205R-020BR 


0.0996 


R (F675W) 


West Limb 


2 X 


700 


1999 Oct 15 


MDM 2.4 m 




0.275 


R 


Whole SNR 


2 X 


720 



^Original Palomar 5 m plate scale: ll'/l per mm. 
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Table 2 

Proper Motion Knot Identification (ID) References 



Knot IDs Prior IDs for Knots Reference 



Outer Knots 



1,4,5,8,9 


1,4,5,8,9 


FBB87, FBG88 


11 


11 


FBG88 


15 


KB91, 15 


vdBK83; FBB87; FBG88 


19 


19 


Fesen 2001 


Knots in NE Jet 






82 


82 


FB91 


99 


115 


vdBK83; FBG88 


116 


116 


KvB76, vdBK83; FBG88 


118,120,121 


118,120,121 


FBG88 


130 


113 


vdBK83 


131 




this paper 


Shell Knots 






93 


41 


KvdB76, vdBK83 


94 


7 


KvdB76, vdBK83 


95 




this paper 


97 


111 


vdBK83 


200 




this paper 


201 




this paper 


302 


102 


KvdB76, vdBK83 


303 




this paper 


304 


101 


KvdB76, vdBK83 


305 


105 


KvdB76, vdBK83 


306 




this paper 


308 




this paper 


401 


4 


vdBD70 


403 


65 


KvdB76 


404 




this paper 


405 




this paper 


41 


59 


KvdB76 


43 


14 


KvdB76, vdBK83 


44 




this paper 


45 


63 


KvdB76 


501 


1 ? 


KvdB76 


981 




this paper 



References: FB91 = Fesen & Becker 1991; FBG88 = Fesen, Becker, 
& Goodrich 1988; FBB87 = Fesen, Becker, & Blair 1987; KvdB76 = 
Kamper & van den Bergh 1976; vdBD70 = van den Bergh & Dodd 
1970; vdBK83 = van den Bergh & Kamper 1983. 
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Table 3 

Knot Positions and Proper Motions 



Knot 




a 


5 




\i6 


fx 




First 


Last 


ID 


(h m s) 


(° ' ") 


(mas yr""*^) 


(mas yr~^) 


(mas) 


(mas yr~^) 


Epoch 


Epoch 


1 


23 


23 12.943 


58 45 57.97 


-341 


-509 


113 


10 


1958, 


.61 


1996.77 


4 


23 


23 47.327 


58 48 26.40 


454 


-70 


172 


9 


1951, 


.83 


1996.76 


5 


23 


23 48.455 


58 48 15.05 


488 


-108 


83 


5 


1951, 


.83 


1999.79 


8 


23 


23 12.202 


58 51 24.06 


-370 


467 


74 


5 


1951, 


.83 


1999.79 


9 


23 


AO Uo . 1 OO 


58 50 46.53 


-458 


362 


69 


4 


1951 


83 


1999.79 


11 


23 


23 11.921 


58 51 16.43 


-391 


454 


143 


15 


1976 


.50 


1999.45 


15 


23 


23 31.727 


58 51 33.51 


97 


500 


77 


4 


1951 


.83 


1999.79 


19 


23 


23 03.424 


58 49 05 79 


— 571 


49 


58 




1976 


.50 


1 qoq 7q 


41 


23 


23 23.993 


58 50 59.85 


-92 


380 


181 


12 


1958 


.61 


1999.79 


43 


23 


23 25.187 


58 50 52.39 


-67 


365 


152 


9 


1951, 


.83 


1999.79 


44 


23 


23 25.367 


58 50 53.47 


-44 


376 


148 


10 


1958, 


.61 


1999.79 


45 


23 


23 24.891 


58 50 55.07 


-71 


369 


150 


9 


1958, 


.61 


1999.79 


82 


23 


23 48.677 


58 50 11.11 


491 


252 


72 


8 


1976 


.50 


1999.79 


93 


23 


23 17.839 


58 50 37.52 


-222 


331 


114 


6 


1951 


.83 


1999.79 


94 


23 


23 18.480 


58 50 39.43 


-225 


327 


97 


5 


1951 


.83 


1999.79 


95 


23 


23 32.801 


58 48 07.52 


110 


-111 


76 


10 


1958 


.61 


1999.79 


97 


23 


23 43.194 


58 48 20.64 


345 


-77 


115 


6 


1958 


.61 


1999.79 


99 


23 


23 52.171 


58 49 53.03 


580 


194 


48 


6 


1958, 


.61 


1999.79 


116 


23 


23 53.316 


58 50 44.52 


604 


358 


78 


6 


1951, 


.83 


1999.79 


118 


23 


23 55.062 


58 50 34.58 


616 


322 


89 


11 


1976 


.50 


1999.79 


120 


23 


23 58.184 


58 50 47.51 


719 


360 


107 


10 


1976 


.50 


1999.79 


121 


23 


23 59.032 


58 51 01.44 


730 


411 


107 


13 


1958, 


.61 


1996.76 


130 


23 


23 49.665 


58 50 08.48 


508 


242 


141 


9 


1958 


.61 


1999.79 


131 


23 


23 53.600 


58 49 57.28 


607 


203 


195 


21 


1976 


.50 


1996.76 


200 


23 


23 41.778 


58 50 03.71 


321 


215 


69 


10 


1976 


.50 


1999.79 


201 


23 


23 44.598 


58 48 33.20 


381 


-4() 


60 


6 


1958 


.61 


1999.79 


302 


23 


23 40.139 


58 50 02.10 


287 


212 


126 


8 


1958 


.61 


1999.79 


303 


23 


23 39.275 


58 50 01.73 


260 


205 


157 


10 


1958, 


.61 


1996.76 


304 


23 


23 39.899 


58 50 07.15 


285 


233 


141 


9 


1958, 


.61 


1999.79 


305 


23 


23 40.732 


58 50 06.15 


307 


225 


162 


10 


1958, 


.61 


1999.79 


306 


23 


23 41.484 


58 50 10.88 


324 


247 


143 


10 


1958, 


.61 


1999.79 


308 


23 


23 41.525 


58 50 08.33 


307 


222 


170 


18 


1976 


.50 


1999.79 


401 


23 


23 22.948 


58 50 04.47 


-108 


219 


124 


8 


1951 


.83 


1999.79 


403 


23 


23 26.438 


58 49 52.17 


-22 


188 


140 


10 


1951 


.83 


1999.79 


404 


23 


23 21.358 


58 50 01.91 


-142 


217 


136 


8 


1958 


.61 


1999.79 


405 


23 


23 20.757 


58 50 05.96 


-172 


232 


158 


11 


1958 


.61 


1999.79 


501 


23 


23 42.343 


58 48 42.61 


329 


-16 


157 


14 


1976, 


.50 


1999.79 


981 


23 


23 47.840 


58 49 43.04 


470 


168 


120 


11 


1951, 


.83 


1976.50 



Note. — Positions are for epoch J2000 and are referred to the ICRS. 
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Table 4 

Age and Center of Expansion Measurements for Cas A 



Reference 


SNR 




To or 




Expansion Center Coordinates 




Region 




SNR Age 




a 


(J2000) 




5(J2000) 


OPTICAL 




















van den Bergh & Dodd 1970 


Brt. Shell 


A.D. 


1667 ±8 


23'' 


23' 


" 2T.16 ±(f.2 


58° 


48' 


47''6 ±3''l 


Kampcr & van den Bcrgh 1976 


Brt. Shell 


A.D. 


1653 ±3 


23'' 


23' 


" 27^76 ±0?1 


58° 


48' 


46''7 ±0''8 




NE ,Ict 


A.D. 


1671 ±3 














van den Bergh & Kamper 1983 


Whole SNR 


A.D. 


1658 ±3 


23'' 


23' 


" 27? 76 ±on 


58° 


48' 


46''4 ±1"0 


Fesen, Becker, & Goodrich 1988 


Outer Knots 


A.D. 


1680 ±15 














Kamper Sc van den Bergh 1991 


Outer Knots 


A.D. 


1671 ±3 














Reed et al. 1995 


Brt. Shell 






23'' 


23' 


" 26! 55 ±0!09 


58° 


49' 


00'.' 7 ±0;'8 


Thorstensen et al. 2001 


Outer Knots 


A.D. 


1671.3 ±0.9 


23'' 


23' 


" 27! 84 ±0!06 


58° 


48' 


49"4 ±0"4 




Shell Knots 


A.D. 


1662.3 ±1.7 


23'' 


23' 


" 27! 71 ±0!09 


58° 


48' 


49'.' 5 ±0;'7 




Whole sample 


A.D. 


1669.1 ±0.8 


23'' 


23' 


" 27! 77 ±0!05 


58° 


48' 


49"4 ±0"4 


RADIO 




















Tuffs 1986 


Brt. Shell 


Age 


= 949 yr 


23'' 


23"" 25?95 ±0!4 


58° 


48' 


48;'4 ±4;'4 


Green 1988 


Brt. Shell 


Age 


~ 400 yr 














Anderson & Rudnick 1995 


Brt. Shell 


Age 


~ 940 yr 


23'' 


23' 


" 26! 05 ±0!2 


58° 


48' 


54:' 3 ±3;'l 




Outer Knots 


Age 


= 550 - 900 yr 














Agiicros Sz Green 1999 


Brt. Shell 


Age 


= 400 - 500 yr 














X-RAY 




















Vink ot al. 1998 


Brt. Shell 


Ago 


= 501 ±15 yr 














Koralesky et al. 1998 


Brt. Shell 


Age 


~ 500 yr 















Table 5 

Derived Center of Expansion versus X-ray Point Source 



Object 


Reference 




Coordinates 










a(J2000) 




(5(J2000) 


Center of Expansion 


this paper 




23m 27S77 ±o!05 


58° 


48' 49'.'4 ±0'.'4 


X-ray Pt. Source 


Chandra ACIS" 


23'' 


23°" 27;86 ±0?13 


58° 


48' 42'.'8 zbl'.'O 


X-ray Pt. Source 


ROSAT HRI'' 


23'^ 


23™ 27!57 ±0!75 


58° 


48' 44'.'0 ±6'.'0 


X-ray Pt. Source 


Einstein HRI* 


23'^ 


23"" 27?83 ±0!50 


58° 


48' 43'.'9 ±4'.'0 


X-ray Pt. Source 


Einstein HRl'' 


23'^ 


23"" 27;89 ±0?50 


58° 


48' 43'.'7 ±4'.'0 



"^Tananbaum (1999); Kaplan, Kulkarni, & Murray (2001); Murray et al. (2001). 
''Prom Pavlov et al. (2000). 
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